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Neuronal morphologies are pivotal for brain functioning: physical overlap between dendrites 
and axons constrain the circuit topology, and the precise shape and composition of 
dendrites determine the integration of inputs to produce an output signal. At the same 
time, morphologies are highly diverse and variant. The variance, presumably, originates 
from neurons developing in a densely packed brain substrate where they interact (e.g., 
repulsion or attraction) with other actors in this substrate. However, when studying neurons 
their context is never part of the analysis and they are treated as if they existed in 
isolation. Here we argue that to fully understand neuronal morphology and its variance 
it is important to consider neurons in relation to each other and to other actors in the 
surrounding brain substrate, i.e., their context. We propose a context-aware computational 
framework, NeuroMaC, in which large numbers of neurons can be grown simultaneously 
according to growth rules expressed in terms of interactions between the developing 
neuron and the surrounding brain substrate. As a proof of principle, we demonstrate that 
by using NeuroMaC we can generate accurate virtual morphologies of distinct classes 
both in isolation and as part of neuronal forests. Accuracy is validated against population 
statistics of experimentally reconstructed morphologies. We show that context-aware 
generation of neurons can explain characteristics of variation. Indeed, plausible variation is 
an inherent property of the morphologies generated by context-aware rules. We speculate 
about the applicability of this framework to investigate morphologies and circuits, to classify 
healthy and pathological morphologies, and to generate large quantities of morphologies 
for large-scale modeling. 
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INTRODUCTION 

Neuronal morphology is important for brain functioning. The 
interplay between dendritic and axonal morphology limits the 
microcircuits (Peters and Payne, 1993), and the shape and com- 
position of dendrites define how inputs are integrated to produce 
outputs (London and Hausser, 2005; Silver, 2010; Torben-Nielsen 
and Stiefel, 2010). As such, it is not surprising that changing mor- 
phological traits and morphological anomalies are implicated in 
neuro-developmental and degenerative diseases (Kaufmann and 
Moser, 2000; Dierssen and Ramakers, 2006). Nevertheless, neu- 
rons come in all shapes and sizes. The diversity is said to express the 
difference between neuron classes while variation represents the 
intra-class differences (Soltesz, 2005). Diversity originates from 
the genetic make-up of neurons (Jan and Jan, 2010; Tavosanis, 
2012). By contrast, the variance can be assumed to originate from 
interactions between the developing neuron and the brain sub- 
strate, its context (McAllister, 2000; Scott and Luo, 2001; Landgraf 
and Evers, 2005; Jan and Jan, 2010; Tavosanis, 2012). Indeed, in 
both axonal (Mortimer et al., 2008) and dendritic (Gao, 2007; Cove 
etal., 2009) development a plethora of microscopic interactions 
have been revealed to influence branching patterns and "guide" 
the direction of growth. Thus, a neuron's context holds the key to 
understanding morphological variance. 

Unfortunately, the context surrounding a neuron has his- 
torically been neglected in the analysis and quantification of 



morphologies. In a highly influential work, Hillman argued that 
dendritic morphologies could be described completely and accu- 
rately by a finite set of morphometric descriptors (Hillman, 1979). 
Thus, the idea was born that careful description of morphometries 
measured from isolated neurons would be sufficient to charac- 
terize neuronal morphology. Later, when digital reconstructions 
became more common practice, this idea inspired the way neurons 
are represented digitally: a pure representation of the morphology 
itself without any information about the context. Currently, a dig- 
ital representation consists of a set of points in three dimensions 
with additional information on how they are linked to each other, 
as is done in the de facto standard SWC format (Cannon et al., 
1998). 

As a consequence, morphometric features used to quantify and 
analyze morphologies (such as the order and degree of points in 
the neuronal tree or neurite lengths) relate to the neuron itself and 
are unable to describe any characteristic of the context. Hence, 
statistical approaches to analyze morphologies and their variance 
that use these morphometric features are bound to fail to describe 
neuronal morphologies correctly as contextual influences includ- 
ing boundaries, capillaries and other neurons, cannot be taken 
into account. Indeed, in earlier work it was shown that the vari- 
ance in morphometric features can be so high that no statistical 
model can be constructed to accurately describe the limited data 
(Torben-Nielsen etal, 2008). 
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An alternative, albeit in practice closely related to the pure sta- 
tistical approach to study neuronal morphologies is the so-called 
"generative approach" (Ascoli etal., 2001; van Pelt and Uylings, 
2002; Stiefel and Sejnowski, 2007; Torben-Nielsen and Cuntz, 
20 14) . In this approach virtual morphologies are generated de novo 
using morphogenetic algorithms. In most cases, these algorithms 
adhere to the ideas proposed by Hillman and sample from statisti- 
cal distributions representing morphometric features to generate 
a morphology (Eberhard et al., 2006; Lindsay et al., 2007; Torben- 
Nielsen and Cuntz, 2014). Clearly, these methods can mimic 
statistical properties of the data set but fail to capture contex- 
tual influences and plausible variation (but see Samsonovich and 
Ascoli, 2003). Notable exceptions exist and target specific char- 
acteristics of the context. Luczak proposed a generative method 
based on diffusion-limited aggregation to illustrate how competi- 
tion over resources and the spatial distribution thereof could shape 
dendritic morphologies (Luczak, 2006). In another work, Cuntz 
and colleagues proposed a generative approach based on high-level 
wiring constraints. By generating multiple virtual morphologies 
in the same volume, competition over resources could be mim- 
icked (Cuntz etal., 2010). In previous work, we demonstrated 
that self-referential contextual cues (e.g., self-avoidance, soma- 
tropism, and membrane stiffness) could be used to explain some 
characteristics of dendritic morphologies (Memelli etal, 2013). 
Recently, CX3D was designed to simulate neuronal development 
based on intrinsic and extrinsic, contextual factors (Zubler et al., 
2013). 

In this work we argue that in order to fully understand neu- 
ronal morphologies we need to break with the view that neurons 
can be treated as independent, isolated entities. Therefore, we 
propose a new approach to study morphologies in which large 
numbers of virtual morphologies are generated simultaneously 
de novo while embedded in a virtual brain substrate, resulting 
in a mechanistic - in contrast to a statistical - description of 
morphologies. In this approach, morphologies are generated by 
repeatedly extending simulated, phenomenological growth cones 
that are guided by interactions with other actors in the brain 
substrate. 

We designed and implemented a prototype of the proposed 
computational framework, NeuroMaC ("Neuronal Morphologies 
and Circuits"). We showcase the functionality of our framework 
related to single neuron morphologies by synthesizing spinal cord 
motor neurons, hippocampal granule cells and cortical layer 5 
(L5) pyramidal neurons. All results are validated against publicly 
available, experimentally reconstructed morphologies. 

MATERIALS AND METHODS 
OUTLINE 

The rationale behind our proposed framework is based on two key 
experimental findings. The first is that the genetic make-up of a 
neuron determines its shape to a large extent. In cell culture exper- 
iments, neurons have a recognizable morphology, albeit one that 
differs from in situ occurrences (Banker and Cowan, 1977; Krieg- 
stein and Dichter, 1984). Second, the genetic make-up of neurons 
also appears to outline a blueprint of neurons in terms of interac- 
tions with the substrate in which they develop. Growth is mainly 
determined by growth-cones that contain filopodia-like structures 



that sense the molecules present in the extracellular matrix. Sen- 
sation of these molecules then influences when a growth cone 
branches or terminates as well as the direction of elongation (Itoh 
etal, 1993; Scott and Luo, 2001; Mortimer etal, 2008; Jan and 
Jan, 2010). 

We extrapolate these key findings to operational concepts in our 
framework that simulates phenomenological growth cones called 
fronts. Broadly speaking, fronts contain growth rules that can be 
expressed in terms of interactions with other agents present in the 
substrate. Interactions are always "local" in the sense that a front is 
able to sample its direct surrounding. As such, fronts are a simple 
metaphor for biological growth cones. 

Figure 1 outlines the concepts underlying NeuroMaC. Based 
on the "local" nature of sensing and sampling of fronts we can 
decompose the simulated brain volume into small sub volumes 
(SVs). Each SV has full knowledge about all contained fronts and 
contextually relevant actors in the substrate, e.g., boundaries and 
other neurons amongst others. All SVs repeatedly extend all active 
fronts contained inside their spanned volume. Because fronts also 
have a physical dimension with a location and a radius, extend- 
ing fronts creates the simulated neurites by creating a frustum 
between the initial position of a front and the new position after 
extension. Details about the construction rules of fronts are pro- 
vided in the next section and for now it suffices to understand 
that - in line with the behavior of growth cones - fronts can 
extend, branch or terminate, and that they can use contextual cues 
to influence these actions. Once the active fronts are extended, 
the SVs perform the crucial step of checking and resolving struc- 
tural overlaps while simultaneously recording locations of putative 
synapses. As a result, generation of morphologies and construc- 
tion of a circuit (without structural overlaps) can be performed in 
one pass. 

NeuroMaC 

We designed and implemented NeuroMaC in accordance to the 
rationale and key concepts outlined above. Here we describe in- 
depth the components of the proposed framework. 

Multi-agent architecture and parallelization 

NeuroMac is designed as a multi-agent system, that is, different 
components of the framework work autonomously and commu- 
nicate with each other through messages. A multi-agent system 
allows straightforward parallelization with the number of com- 
puting cores to ensure scalability. NeuroMaC has two agent types: 
one central administration agent and multiple SV-agents. 

The administration agent performs all internal housekeeping. 
It reads a configuration file (Table 1) that defines the simula- 
tion and system specific settings. Subsequently the administration 
agent decomposes the brain substrate into smaller SVs and ini- 
tializes the SVs. During initialization each SV is assigned a space 
it controls together with all environmental details required for 
the fronts to develop. The administration agent maintains a cen- 
tral clock to synchronize updating of fronts in each SV. A clock 
ensures that irrelevant issues such as execution time on the com- 
puting resource do not bias simulated growth. In case an updated 
front moves outside of the space covered by a particular SV, the 
administration agent brokers the migration of that front to the 
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FIGURE 1 | Schematic of the proposed context-aware framework, 
NeuroMaC, to generate virtual morphologies. (A) The simulated brain 
substrate is decomposed into small sub volumes (SVs). Sub volumes keep 
track of all neurites and other relevant actors inside their spanning volume. 
(B) Algorithm performed by each sub volume during one simulated, 
centrally controlled time step. (C) Fronts are implemented as cellular 
automaton-like machines and conceptually related to growth cones in that 



they update their location based on the local context. Full lines: neurites 
(black and gray: existing; green: newly added). Circles represent active 
(filled) or inactive (open) fronts. Dashed lines represent the contextual cues 
influencing the direction of growth of an active front to be extended 
(indicated by a red circle). Here the contextual cues are defined by an 
inertial forward-directed influence, another neurite, and a gradient in the 
substrate. 



appropriate new host SV. All updates inside an SV are communi- 
cated to the central agent, which compiles a centralized output file 
containing all neuronal morphologies. 

The SV agents perform the same behavior in parallel. The 
number of these agents can scale with the number of available 
computing nodes; more nodes results in smaller decomposed 
volumes and faster run times. Conceptually, SVs represent the 
direct neighborhood surrounding a developing growth cone. Dis- 
tal parts of the brain substrate are of no concern to a growth 
cone as all contextual cues are sensed in the direct vicinity. SVs 
contain all local information about the substrate itself, e.g., bound- 
aries, laminar structure, same and other neuron structures, etc. 
Diffusible molecules in the extracellular space can promote long- 
distance interactions and while we do not simulate diffusion 
explicitly, the effect of contextual cues can propagate from SV 
to SV so that these are also locally available for growth cones. 
Any cue not on the hosting SV or on one of the direct neighbors 
is summarized (averaged) and only this information is revealed 
to active fronts. This measure is valid because it is irrelevant 
for an active front to know the exact locations of very distant 
cues. 

During each general time step SVs execute the algorithm listed 
in Figure IB. However, just before the algorithm is executed, each 
SV communicates with its neighbors to query their contained vol- 
ume. This is needed because, if an active front is close to an SV 
boundary (e.g., close enough that it might interact with a neurite 
contained in a neighboring SV), it also has to sense the neigh- 
boring substrate. During the main algorithm, SVs call each active 
front inside their volume, in randomized order, to compute its 



next location (see next section). Once the SV receives the updated 
front, it performs several checks. First, it checks if the new loca- 
tion of the front is still inside the volume it spans. If not, the 
front is migrated to another SV. Otherwise, the SV checks whether 
the new front physically overlaps with existing fronts and neu- 
rites. Overlap is tested between two fronts and their associated 
line segments. That is the line segment between a front and its 
parent. If the minimal distance between two such line segments 
is smaller than the sum of the radii of both associated fronts we 
consider this to be an overlap. Unless the radius of a front is dras- 
tically smaller than that of its parent front, this method yields 
adequate results. When a potential overlap is detected, the SVs 
will try to resolve it by randomly perturbing the front's location. 
If the conflict cannot be resolved in a predetermined number of 
attempts, the front is terminated at its previous position. When all 
active fronts are updated and validated, the corresponding newly 
formed neurites are communicated to the administration agent. 
Putative synapse locations are computed in the same way (and at 
the same time) as the structural overlaps with the difference that 
a maximally allowed distance is set by the user that reflects the 
pre-synaptic bouton and post-synaptic spine size. Although rudi- 
mentary, this method yields a list of putative synapse locations 
that can be pruned in a post-processing step (Hill et al., 2012), but 
also see van Pelt et al. (2010). 

Growth cones as cellular automata 

In NeuroMaC fronts are phenomenological implementations 
resembling biological growth cones. An active front is a front that 
is still developing; an inactive front becomes continuation point, 
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Table 1 I Exemplar configuration file used in NeuroMaC. 



[ system] 

# framework related settings 
seed — 2 

proxy_sub_port = 5599 
proxy_pub_port = 55 60 
pull_port = 55002 
time_out = 10000 

# simulation related settings 
no_cycles — 105 

out_db — models/L5_pyramid/f orest_Z8 . db 
synapse_distance — 5 

# attempts to resolve overlap-conflicts 
avoidance_attempts — 2 

[substrate] 

# settings about the simulated brain substrate 
dim_xyz = [6000.0,1800.0,1410] 

# volume decomposition into xa x ya x za SVs 
xa=6 

ya=6 
za— 1 

# laminar structure 

virtual_LAYER = {6: [[0,0,0], [2000, 2000, 471]], \ 

5 : [ [0, 0,471], [2000, 2000, 826]], \ 
4:[[0, 0,826], [2000, 2000, 1090]], \ 
3: [ [0, 0, 1090], [2000, 2000, 11 92 ]],\ 
2:[[0, 0,11 92], [2000, 2000, 1311]], \ 
1 : [ [0, 0,1311], [2000,2000,1406]]} 

#pia as boundary 

pia — models /L5_pyramid/pia_f orest . pkl 
[cell_type_l] 

# settings related to the growth rules 
no_seeds — 100 

algorithm — Full_detailed 

location = [ [250 , 250 , 800 ],[ 5750 , 1550 , 1 180 ] ] 
soma_radius — 10 



The configuration follows the Python ConfigParser structure. Parameters are 
pooled in several sections and parameter values can take the form of executable 
Python statements. A description is in the main text. 

branching point or a terminal tip. As such, neurites are represented 
by frusta connecting subsequent fronts (Figure 1C; Cannon etal., 
1998;Ascoli etal, 2007). 

Fronts have a dual identity. On the one hand they are physi- 
cal structures with a location and radius in space. On the other 
hand, a front is a cellular automaton-like machine that contains 
its own growth rules describing how and when it should extend, 
branch or terminate (see Table 2 for an example). When an active 
front is not terminating, it either produces one or two new fronts; 
the old front becomes inactive and the newly formed front(s) 
become(s) active fronts. The location of the new front is com- 
puted in accordance to a front's construction rules and locally 
available information. Information can be everything that is con- 
tained in the SV. For instance, homotypic (Grueber etal., 2005; 
Marks and Burke, 2007; Memelli et al., 2013) and same-type (Scott 
and Luo, 2001; Jan etal., 2003) cues can be used, or, the transient 



laminar information through which a front might travel (Hevner 
etal., 2003; Chen etal., 2005). The aforementioned cues have a 
direct biophysical interpretation, but also more phenomenologi- 
cal cues such as (directional) information related to a boundary 
can be used in our framework. A biological counterpart thereof 
could be envisioned to be Reelin secreted by Cajal-Retzius cells 
(Frotscher, 1998; Marin-Padilla, 1998). Construction rules define 
how the front interacts with these other inhabitants of the SV: no 
interaction, repulsion or attraction. Hence, the context is used 
as a guidance cue (Figure 1C). The influence of these cues can 
be distance-dependent mimicking gradients of secreted molecules 
(Mortimer etal., 2008). In addition, fronts can also modify the 
substrate by secreting entities: phenomenological representations 
of secretion molecules that can in turn be used as a guidance cue 
(Hentschel and van Ooyen, 1999). 

IMPLEMENTATION 

We implemented a prototype in Python and use ZeroMQ 
(Hintjens, 2013) to send messages between the components 
because it has the ability to buffer large messages and operate 
asynchronously. The algorithm underlying the behavior of an 
active front is a Python script and is the only part that has to 
be implemented by an end-user. This prototype is available on 
https://groups.oist.jp/cnu/neuromac. 

Combined, the eminent features of NeuroMaC are: (1) 
Context-aware generation of virtual morphologies that will not 
overlap with one another in space; (2) The ability to detect and 
record synapses on the fly; and (3) Straightforward scalability and 
parallelization to generate large numbers of morphologies at the 
same time. 

RESULTS 

In order to validate the proposed framework we generated sets of 
virtual neuronal morphologies and compared them to the statis- 
tics of experimentally reconstructed morphologies. We validate 
NeuroMaC by demonstrating that we can (1) generate mor- 
phologies in isolation as current state of the art approaches do, 
(2) populate a space by generating a forest of non-overlapping 
and interacting hippocampal granule cells, and (3) generate fully 
context-aware morphologies that interact with the environment 
(L5 pyramidal neurons in a laminar architecture). We selected 
these neuron types because motor neurons and hippocampal 
granule neurons are often used in algorithmic generation; pyra- 
midal neurons were chosen because their higher morphological 
complexity and assumed context-dependence. The experimentally 
reconstructed neurons were downloaded from NeuroMorpho.org 
(Ascoli etal., 2007). We took two motor neuron archives, the 
Burke archive (N = 6, Cullheim etal., 1987) and the Fyffe 
archive (N = 8, Alvarez etal, 1998). The granule neurons come 
from the Lee archive (N = 25, Carim-Todd etal., 2009). Pyra- 
midal neurons are layer 5, secondary motor cortex neurons 
and come from the Kawaguchi archive (N = 10, Hirai etal., 
2012). 

MOTOR NEURONS IN ISOLATION 

Motor neurons have a relatively straightforward morphology 
that, from the point of view of an external observer, is fairly 
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Table 2 | Complete Python code used to implement the growth rules underlying the generated motor neurons (illustrated in Figure 2). 



from growth_procs import unit_sample_on_sphere, \ 
direction_to, \ 
gradient_to, \ 
normalize_length, \ 
get_entity, \ 
get_eigen_entity, \ 
prepare_next_f ront 

L_NORM — 40 # fixed-size elongations 

def extend_front (front, seed, constellation) : 

if front . order—— 0 : # this is the soma, create the stems 
new_f ronts — [ ] 

for i in range ( np . random . randint ( 8 , 1 7 )) : 
rnd_dir — unit_sample_on_sphere ( ) 

new__pos — front.xyz + normal! ze_length (rnd_dir, L_NORM) 
new_f ront — prepare_next_f ront (front, new_pos, \ 

set_radius— 8 . 0, add_order— True ) 
new_f ront . swc_type— 2 
new_f ronts . append (new_f ront ) 
return new_fronts 
else : 

# Follow a simple branching rule in all other cases 
bif_prob — 0.6 / ( front . order *2 . 5 ) 
if front . order > 5 : 
bif_prob =0.03 

if np . random . random ( ) > bif_prob: # continue a front 

# random component 

rnd_dir = unit_sample_on_sphere ( ) 

# unit vector of current heading 

heading=normali ze_length ( front . xyz - front. parent. xyz , 1 ) 

# soma-tropism, sample direction away from the soma 
soma_dir = -1.0 * normali ze_length ( direct ion_to (front , \ 

[ front . soma_pos ] , what = "nearest") ,0.4) 

# combine all infliences on the new direction of growth 
new_dir = heading, 1 . 0+ soma_dir + rnd_dir 

new__pos = front.xyz + normali ze_length (new_dir, L_NORM) 
new_front — prepare_next_f ront (front, new_pos, \ 

radius_f actor = 0 . 9 , add_order = False) 

if np . random . random ( ) < 0.06 and front . path_length > = 600: 
return [] 

return [ new_f ront ] 
else: # branch a front, generate two child fronts 

new_f ronts = [] 

for i in range (2) : 

rnd_dir = unit_sample_on_sphere ( ) 

heading = front.xyz - f ront . parent . xyz 

new_dir = normalize_length (heading, 1 . 5) +rnd_dir 

new_pos = front . xy z+normalize_length (new_dir , L_N0RM) 

new_front = prepare_next_f ront (front, new_pos, \ 

radius_f actor = 0 . 7 , add_order = True) 
new_f ronts . append (new_f ront ) 

return new_fronts 



NeuroMaC contains auxiliary function to build fronts and sample the context; these functions are first imported. The main function "extend_front" is called by the sub 
volume and contains the actual growth rules. In this example, a single contextual cue, soma-tropism, is used. 



context-independent (Figures 2A-C). We devised a purely phe- 
nomenological growth rule to mimic the final morphology con- 
sisting of two sub rules: one rule for the initial front (=the 
soma) and one rule for all other fronts. The full Python 
code of the growth rule is listed in Table 2. At the soma 
("front.order == 0"), multiple stems are created in random 
directions around the soma. Once the stems are created fronts 



can bifurcate with a probability inversely proportional to the 
branching order, terminate with a small probability or extend 
otherwise. When a front grows outside the assigned substrate 
space it is terminated. Current heading, repulsion by the soma 
and a random component set the direction of a bifurcating or 
extending front. Typical resultant virtual morphologies are listed 
in Figures 2D-F. 
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FIGURE 2 | Validation of generated alpha motor neurons. (A-C) 

Exemplar experimentally reconstructed spinal cord alpha motor neurons 
|A,B from the Fyffe archive (Alvarez etal., 1998), C from the Burke archive 
(Cullheim etal., 1987)]. (D-F) Virtual morphologies generated by 
NeuroMaC. (G-l) Quantitative comparison. Population morphometries are 
shown for the Burke ("Burke") and Fyffe ("Fyffe") archives and for the 



generated morphologies ("Syn"). (G) Euclidean distance between the 
soma and each terminal point in all morphologies. (H) Topological order of 
each branching point in all morphologies. (I) Occurrence of branching 
points in each morphology as a function of Euclidean distance (i.e., 
Sholl-intersections, see main text). See Table 3 for detailed statistics of 
these (and other) morphometries. 



Visual inspection shows high resemblance between the exem- 
plar and generated motor neuron morphologies. We then checked 
the global morphometric, namely the Euclidean distance between 
the soma and terminal tips (Figure 2G) and the two-dimensional 
local metrics "order" that expresses the occurrences of branch- 
ing points as a function of branching order (Figure 2H), and, 
"Sholl-like," a quick implementation of the Sholl metric that mea- 
sures branch points as a function of Euclidean distance from 
the soma (Figure 21). Trends contained in the experimentally 
reconstructed neurons (labeled "Burke" and "Fyffe") are repli- 
cated by the generated neurons (labeled "Syn"). We quantify 
the distribution by the median (M) and median absolute devi- 
ation (MAD) because the shape of the resultant distribution 
of the measures is unknown a priori and do not necessar- 
ily follow a normal distribution. Spread of the distribution is 
quantified with the interquartile range (IQR). Quantification 
is listed in Table 3. From the quantification we can see that 
there is a fair difference between the exemplar archives and 
that the generated neurons fit well between the values of the 
exemplars. 

Both visual inspection and the quantitative measures show a 
good correspondence between the experimentally reconstructed 
and generated morphologies. These results are on par with the pre- 
viously published results (Memelli et al., 2013), and hence we can 



conclude that by using NeuroMaC we can create sets of neurons 
generated in isolation. 

A FOREST OF HIPPOCAMPAL GRANULE NEURONS 

Next we set to generate granule cells, both in isolation and in a 
"forest" setting, that is, many neurons packed in one volume with 
all neurons being generated simultaneously. Three experimentally 
reconstructed exemplar morphologies are shown in Figures 3A-C. 
We devised a straightforward construction rule in a similar vein to 
the rule used for the virtual motor neurons. Once the soma and two 
initial branches are created, branching occurs with a probability 
that decreases with the centrifugal order of the front. The direction 
of growth is determined by repulsion away from same-neuron 
dendrites, the current heading of a dendrite, and the direction 
towards the superficial part of the substrate, which in this case is 
the superficial part of the dentate gyrus. A random component 
is added to all growth directions as well. We generated two sets 
of virtual morphologies, namely a set in which each neuron was 
generated in isolation (N = 25, Figures 3D-F are representative 
examples) and one set in which 100 morphologies were generated 
simultaneously in a (Figure 3G). The growth instructions were 
kept identical in both sets. The simulated volume, however, was 
increased 20-fold in the forest setting (i.e., 1300 |x x 300 u, x 225 [L, 
with 225 |x being a plausible depth of the dentate gyrus). Note that 
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Table 3 | Quantification of generated and experimentally 
reconstructed alpha motor neurons. 







Synthetic 


Burke 


Fyffe 


# branch points 


M 


125 


161 


55 




MAD 


11 


12.5 


20 




IQR 


41 


23.7 


82 


Euclidean D 


M 


750 


926 


616 




MAD 


171 


184 


277 




IQR 


352 


372 


380 


Max order 


M 


7 


9 


7 




MAD 


0 


0 


1 




IQR 


1 


0 


2.5 


Order 


M 


3 


4 


4 




MAD 


1 


1 


1 




IQR 


2 


2 


3 


Sholl-like 


M 


600 


575 


445 




MAD 


240 


171 


213 




IQR 


480 


350 


435 


Total length 


M 


69,674 


105,373 


28,876 




MAD 


14,041 


8730 


14,691 




IQR 


24,220 


14,169 


41,363 



Distribution of observed morphometries are given by the median (M), median 
absolute deviation (MAD) and inter-quartile range (IQR). Values shown for the 
generated ("synthetic") morphologies and the morphologies originating from the 
Burke and Fyffe archives (see main text). 

in the "forest" setting, developing morphologies interact indirectly 
with each other through overlap-prevention. 

Visually the generated morphologies bear strong resemblance 
to the exemplar ones. We then measured the Euclidean distance 
between some and terminal tips and the maximum order in a tree 
(Figures 3H,I), as well as the two-dimensional "Order" and "Sholl- 
like" metric (Figures 3J,K) for the set of exemplar morphologies 
("Lee") and the sets of morphologies generated in isolation ("Syn") 
and in a forest setting ("Forest"). To avoid biases introduced to an 
unequal number of samples, we randomly picked 25 morpholo- 
gies from the forest and computed the appropriate features from 
this subset. The histograms indicate similar trend in the data of all 
data sets. Quantification of all measured morphometries is pro- 
vided in Table 4. It is interesting to note that the variance in the 
morphologies generated in a forest setting is higher. This obser- 
vation results from the fact that all these neurons are generated 
simultaneously. As a result, some branches would overlap with 
each other. Overlaps are detected and an attempt is undertaken 
to resolve the overlap. However, if no quick resolution is found, 
the branch is terminated. In the forest setting, the somata are 
close to each other and some conflicts in the proximal branches 
could not be resolved and caused very small Euclidean length 
and low maximal order in rare cases (Figures 3H,I, left-most 
red bars). The two dimensional metric indicate a good match 
in the topological and geometrical distribution of branch points 
(Figures 3J,K). 



Even though the neurons in the forest setting were densely 
packed (Figure 3G) no overlaps occurred as neurite locations 
were either corrected or terminated during the validity checks 
performed by the SVs. Therefore, we conclude that with Neu- 
roMaC we can generate forests of non-overlapping, plausible 
morphologies. 

CONTEXT-AWARE L5 PYRAMIDAL NEURONS 

As a final demonstration of the capabilities of NeuroMaC, we 
generated context-dependent layer 5 pyramidal neuron mor- 
phologies. Three exemplar morphologies are shown in Figure 4A. 
By visually examining these morphologies, we can observe some 
morphological traits such as a difference in "height" but these traits 
are hard to relate to their context. However, from canonical circuit 
information, we know that the somas are located in layer 5, that 
their basal dendrites remain mainly in L5 and may extend a bit 
into L4, that their apical dendrite extends to the superficial parts 
and ends close to the pia (in LI) after branching extensively in 
layers L3-L1, and, that oblique dendrites sprout from the apical 
trunk in L4. The remarkable difference in "height" of the apical 
tree, is a clear signature of this context dependence as more super- 
ficially located pyramidal cells cannot extend as far as more deeply 
positioned ones. 

We designed construction rules that take these canonical, 
contextual traits based on laminar structure into account. A 
truncated code snippet is listed in Table 5 to indicate particu- 
lar context-dependent growth rules. Note that the growth rules 
are different for basal and apical dendrites, and a further divi- 
sion of the apical growth rules into rules for L5/L4, oblique 
dendrites, and the dendrites in L3/L2/L1. At the soma, we gen- 
erate an appropriate number of basal stems and one apical stem. 
The basal dendrite branches with a probability inverse propor- 
tional with the centrifugal order; at orders higher than 6 no 
branching is allowed. Termination of a basal branch occurs with 
a small probability or when a branch grows outside the lim- 
iting volume. Direction of growth is again influenced by the 
heading and same-neuron repulsion and an additional random 
factor. The apical branch is contextually aware and the con- 
struction rules change depending on the layer it is in (Table 5, 
"extend_apical_front"). Layer-dependent behavior is biologically 
feasible because in cortex some transcription factors are exclu- 
sively expressed in layer specific neurons (Hevner et al., 2003; Chen 
etal, 2005). In layers 5 and 4, oblique dendrites can sprout and 
grow away from their initial branch point at the apical trunk. In 
subsequent layers (3, 2, and 1) neurons can branch with layer 
specific probabilities as long as a maximum increase in order 
has not occurred yet in one layer. Same-neuron repulsion, cur- 
rent heading, a distance-dependent attraction to the pia, and a 
random component determine the direction of growth in the 
superficial layers 3—1. Apical neurites can terminate as soon as 
they reach layer 3 (and later 2 and 1) with a small probability. 
All apical neurites are terminated if the pia is closer than 35 |x 
away. 

Two sets of morphologies are generated; again one with neu- 
rons in isolation (N = 10 to match the sample size in the 
Kawaguchi archive) and one with 100 simultaneously generated 
morphologies in a forest setting. The volume in the "forest" setting 
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FIGURE 3 | Validation of generated hippocampal granule cells. (A-C) 

Experimentally reconstructed granule cells (from the Lee archive; Carim-Todd 
etal., 2009). (D-F) Virtual morphologies generated by NeuroMaC. (G) Forest 
of 100 simultaneously generated, non-overlapping granule cells. (H-K) 
Quantitative comparison. Population morphometries are shown for the Lee 
archive ("Lee"), synthetic neurons generated in isolation ("Syn") and as part 



of a forest ("Forest"). (H) Euclidean distance between all terminal tips and 
the soma. (I) Maximum topological order in the individual morphologies. (J) 
Topological order of each branching point in all morphologies. (K) Occurrence 
of branching points in each morphology as a function of Euclidean distance 
(i.e., Sholl-intersections). See Table 4 for a detailed quantification of these 
(and other) morphometries. 



was a rectangle of size 6000 u. x 1800 [i x 1400 [i, where 1400 u. is 
the estimated depth of L5 in the exemplar data. All morphologies 
from the former set are plotted in Figure 4B along with the canon- 
ical virtual laminar architecture in which they grew (blue line: pia, 
red dashed lines: layer boundaries. Layer 1 is at the top and layer 5 
at the bottom; layer 6 is not shown). The forest from the latter set 
is plotted in Figure 4C. 

Visually, the generated neurons clearly exhibit the morpholog- 
ical traits summarized above. Furthermore we compared the total 
number of branch points (Figure 4D), the Euclidean distance to 
the terminal tips (Figure 4E) and the total length (Figure 4F). A 
quantification of all measured morphometries is listed in Table 6. 
The basal and apical dendrites are treated separately in these 
measures. The basal trees show great correspondence with the 



exemplar morphologies in terms of the Euclidean distance to 
the terminal tips and the total length of the dendritic trees. The 
number of branch points in the generated neurons is markedly 
higher than in the exemplar ones; a range of [19,39] for the 
Kawaguchi archive and [20,52] and [19,53] for the generated 
neurons in isolation and forest setting, respectively. Given a 
correct match with the total length and the Euclidean distance 
to the tips, we speculate that the simple branching and ter- 
mination rules are not sufficient for the basal trees, although 
the low number of branch point can also result from incom- 
plete reconstructions (Anwar etal., 2009, but also see Section 
"Discussion"). 

Considering the apical trees, we observe a mismatch in the 
Euclidean distances and the total length between the exemplar 
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Table 4 | Quantitative description of experimentally reconstructed 
hippocampal granule neurons and their generated counterparts. 







Lee 


"Isolation" 


"Forest' 


# branch points 


M 


13 


12 


13 




MAD 


1 


4.5 


5 




IQR 


2 


8.25 


9 


Euclidean D 


M 


207 


199 


197 




MAD 


12.7 


11.7 


16.7 




IQR 


26.2 


41.7 


68 


Max order 


M 


5 


5 


5 




MAD 


0 


1 


1 




IQR 


1 


1 


1 


Order 


M 


3 


3 


3 




MAD 


1 


1 


1 




IQR 


1 


1 


1 


Sholl-like 


M 


77 


79 


75 




MAD 


38.5 


40.5 


40.1 




IQR 


84 


85 


80 


Total length 


M 


2255 


1846 


1590 




MAD 


258 


437 


505 




IQR 


391 


1001 


922 



Generated morphologies can be generated in isolation or in a forest setting. Pre- 
sentation as in Table 3 Values shown for the morphologies from the Lee archive 
and the morphologies generated in "Isolation" and in the "Forest" settings (see 
main text). 

and the generated morphologies. We attribute both to a dif- 
ference in the oblique dendrites. As seen in Figure 4E (left 
panel, "Kawaguchi"), there is a peak of terminals in the apical 
dendrite that terminate close to the soma. While the gener- 
ated data also displays a second peak due to terminals of the 
oblique dendrites, this peak is less pronounced and shifted to 
greater Euclidean distances. We speculate that in the exemplar 
dendrites, more oblique dendrites sprouted more proximally than 
in our model. Given a major thalamic synaptic pathway in cor- 
tex projecting to layer 4 and synapsing onto oblique dendrites 
(Meyer etal, 2010; Oberlaender etal., 2012), it is not unreason- 
able to think the oblique dendrites mainly sprout in layer 4 as 
in our model. But, as said, an SWC file does not contain any 
contextual information so the true dimensions of the laminar 
architecture of the animals from which the neurons were recon- 
structed remain a guess. Moreover, we consider the ability of 
NeuroMaC to construct context-dependent dendrites a quality, 
even if no context-dependent information related to the exemplar 
morphologies was directly available. The fact that the apical trees 
generated by NeuroMaC all reach the LI - and not further - are a 
great illustration of this context-dependence. 

Our results indicate a clear and valid context-dependence, 
which is similar to the morphological traits in the exemplar 
data. Therefore, we can conclude that the generated morpholo- 
gies exhibit context-dependent morphological traits that match to 
the traits discovered in the exemplar data. 



DISCUSSION 

We started this work with the observation that there is a large 
discrepancy between the way neuronal morphologies are studied 
(in isolation) and the way they develop and take their shape (in 
interaction with a dense surrounding substrate). From experi- 
mental studies it appears that the surrounding brain substrate, the 
context of all neurons, plays a pivotal role in shaping the morphol- 
ogy and resultant brain circuits. To overcome this discrepancy, we 
proposed a new computational framework, NeuroMaC, to study 
how neuronal morphologies emerge from interactions with other 
actors in the brain substrate. 

We opted for a phenomenological framework for the sake of 
conceptual simplicity and to curb computational costs. Construc- 
tion rules are conceptually related to the genetic make-up of a 
neuron and express how a neuron has to grow in terms of repul- 
sive or attractive interactions with the surrounding substrate. 
A phenomenological framework helps to reduce the computa- 
tional resources in contrast to biologically and physically detailed 
ones. Moreover, the design of NeuroMaC as a multi-agent system 
ensures scalability with the number of available processors. As a 
consequence of the design choices, NeuroMaC can be used to gen- 
erate large numbers of interacting morphologies simultaneously. 
This feature is unrivaled. CX3D, an existing computational tool 
aims to simulate the whole of cortical development, from migra- 
tion over polarization and differentiation to dendrite and axon 
formation. However, the main version is serial (i.e., not parallel) 
which limits its applicability to generate multiple full morpholo- 
gies at the same time. NETMORPH, a tool capable of generating 
large cortical networks (Koene etal, 2009) adopts a strategy in 
which a volume is populated by adding neurons that are generated 
in isolation. The topology of neurons is based on a mechanis- 
tic growth rule but the geometry assigned to embed the topology 
in space is statistically sampled from exemplar data. Hence, in 
NETMORPH all neurons are independent and not based on any 
contextual cues (van Ooyen etal., 2014). Although it has to be 
noted that exemplar data contains morphologies that are shaped 
through contextual interactions and, therefore, if a model suc- 
ceeds in reproducing morphological traits it implicitly captures 
some of these interactions. Historically, ArborVitae (Senft and 
Ascoli, 1999) was proposed to generate large networks of neu- 
rons simultaneously and with some phenomenological interaction 
based on resource competition. While promising initial results 
were generated, this tool is no longer in development. Hence, Neu- 
roMaC is currently the only computational framework to study 
explicitly how neurons grow together while interacting with the 
environment. 

We demonstrated that by using NeuroMaC we can generate 
plausible neuronal morphologies with construction rules based 
on local interactions, which inhabit the same simulated substrate 
and have no physical overlaps. In the current work, construction 
rules underlying the growth of morphologies are a crude approx- 
imation of the hypothesized growth rules used by neurons. The 
aim of this work was not so much the generation of the most "real- 
istic" morphologies or morphological traits but rather showcasing 
the power and usability of our new framework. As such, we illus- 
trated that construction rules expressed in terms of repellants and 
attractors are a useful metaphor to study morphologies. 
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FIGURE 4 1 Validation of generated layer 5 pyramidal neurons. 

(A) Experimentally reconstructed layer 5 pyramidal neurons (from the 
Kawaguchi archive). (B) Virtual morphologies generated by NeuroMaC. 
Simulated laminar structure (L1-L5, from top to bottom) indicated by 
dashed lines; blue line represents the pia. (C) Forest of 100 
simultaneously generated, non-overlapping pyramidal neurons. (D-F) 
Quantitative comparison. Population morphometries are shown for the 



Kawaguchi archive, synthetic neurons generated in isolation ("Syn") 
and as part of a forest ("Forest"). Statistics are given for basal 
(left panels) and apical (right panels) trees separately. Shown are 
total number of branching points (D), Euclidean distance between 
terminal tips and the soma (E) and the total length of the 
dendrites (F). Detailed statistics of these (and other) morphometries 
in Table 6. 



NeuroMaC can be used in any desired way on the contin- 
uum between small and large spatial scales and their associated 
level of biological detail. At one end of this continuum it can 
be used to study the effects of detailed, biologically plausible 
construction rules. This way, studies can be conducted investi- 
gating how particular construction rules representing biophysical 
processes influence morphological traits. On the other end of 
the continuum, one could opt to use less detailed rules to 
generate full morphologies and, because putative synapse loca- 
tions are recorded as well, the resultant circuits. Of course, 
highly detailed construction rules can also be used (at little 
extra computational cost) to generate full circuits and any "inter- 
mediate" level of detail can be implemented as well. However, 
while it is possible to compute the propagation of micro- 
scopic rules to the meso-scale circuit, it can be a tedious 



task to analyze the whole circuit at large for traces of the 
underlying microscopic interactions. Another noteworthy fea- 
ture of NeuroMaC is that it supports a mixed-methodology 
with respect to the growth rules. That is, existing context- 
independent neurogenetic algorithms can be implemented in 
a straightforward fashion so that they can be used as growth 
rules. As such, a simulated brain substrate could be populated 
by morphologies grown in accordance to different methodolo- 
gies. 

One important observation is that our virtual morphologies 
generated in a forest setting exhibit a larger variance than present 
in the exemplar data (Figures 3IJ and 4D,F). This effect is smaller 
but still present in the neurons generated in isolation. We turn to 
the data sets of experimentally reconstructed neurons to explore 
the issue of variance. 
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Table 5 | Code snippet illustrating the growth rules to generate layer 5 pyramidal neurons. 



def extend_f ront (front, seed, constellation) : 
if front . order — — 0 : 
new_fronts — [] 

apical_front — create_apical_branch (front, constellation) 
basal_f ronts — create_basal_branches (front, constellation) 
new_f ronts . append (apical_f ront ) 
new_f ront s . extend (basal_f ront s ) 
return new_f ront s 
elif front . swc_type — — 3 : 

i f front . update_cycle < — np . random . randint (35, 47) : 

return extend_basal_f ront (front, constellation) 
else : 

return [ ] 

else 

return extend_apical_f ront (front, constellation) 
def create_apical_branch (front, constellation) : 

# create one branch in direction of the pia 
pi a — get_entity ( "pia" , constellation) 

dir_t o_pia — direct ion_to ( front , pia, what = "nearest " ) 
new_dir = normal i ze_length ( dir_to_pia, 3.0) 

new_pos = front . xyz + normalize_length (new_dir , APICAL_NORM) 
new_f ront — prepare_next_f ront (front , new__pos , \ 
set_radius = 1.0, add_order = True ) 
return new_f ront 

def create_basal_branches (front, constellation) : 
for i in range ( np . random .randint ( 5 , 11) ) : 

# construct a number of basal branches 

def extend_basal_f ront (front, constellation) : 

# branch, continue or terminate 

def extend_apical_f ront (front, constellation) : 

# terminate branches too close to the pia 

# sprout oblique dendrites in L5 or L4 

if (front. layer — — 4 or front . layer — — 5 ) and not front. oblique : 

# special rule for oblique dendrites 
if front . oblique : 

# compute next location and return 

# cont inue or terminate oblique branch? 

# layer specific rules for front s in different layers 
if front. layer > — 3: 

# branch, continue or terminate 

if front . layer — — 2 : 

# branch, continue or terminate 
if front. layer — — 1: 

# branch, continue or terminate 



The code is incomplete and merely for the purpose to illustrate some of the context-dependent cues such as growth direction to the pia and layer specificity (for the 
apical tree). 



We can start by assuming that the data is a good represen- 
tative of all neurons. In that case, our data exhibits too much 
variation. Here the explanation would be that the used branching 
rules are too simple and that branch probability and termina- 
tion are also dependent on both intrinsic and extrinsic signals. 
Intrinsic signals could be mediated through the production and 
transport of actin filaments that are required for scaffolding the 
neuronal membrane (Graham and van Ooyen, 2004). A detailed, 
mechanistic rule based on these intrinsic properties has been pro- 
posed (van Pelt and Verwer, 1986; van Pelt and Schierwagen, 



2004) and could be used in our framework. Extrinsic signals 
are inherently context-dependent. Experimental work has demon- 
strated that the presence of specific molecules in the extra- cellular 
space influence branching and termination properties (Itoh etal., 
1993; Dimitrova etal., 2008). While we did not address biologi- 
cally plausible termination and branching conditions, we did use 
the contextual laminar architecture as a cue to set layer specific 
branching probabilities, and fronts in close proximity to the pia 
were terminated. Another way of restricting virtual morphologies 
is by generating them inside a limited space as applied here to 
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Table 6 | Quantitative description of experimentally reconstructed L5 
pyramidal neurons and their generated counterparts. 



Kawaguchi "Isolation" "Forest" 



TT Dl al IL.I 1 


3 ni^a 1 
apiud 1 


M 


46 


40 5 


33 


nnint^ 

1 ! 1 Lo 




MAD 


3 


10 


5 






IQR 


4.7 


13.7 


10 




basal 


M 


31 


30 


32 






MAD 


1.5 


6 


5.5 






IQR 


2.5 


14.2 


10 


Euclidean D 


apical 


M 


407 


561 


552 






MAD 


198 


116 


98.3 






IQR 


409 


265 


237 




basal 


M 


125 


138 


137 






MAD 


28.5 


35.7 


38.8 






IQR 


62.1 


76.7 


85.9 


Max order 


apical 


M 


17.5 


19 


18 






MAD 


2.5 


1 


2 






IQR 


4.7 


2.5 


4.7 




basal 


M 


5 


5 


5 






MAD 


1 


0 


0 






IQR 


1.7 


0 


0 


Order 


apical 


M 


10 


16 


15 






MAD 


5 


3 


3 






IQR 


10 


8 


8 




basal 


M 


2 


3 


3 






MAD 


1 


1 


1 






IQR 


1 


2 


2 


Sholl-like 


apical 


M 


345 


521 


554 






MAD 


224 


171 


131 






IQR 


448 


378 


322 




basal 


M 


41 


70 


70 






MAD 


13.3 


35 


42 






IQR 


28.3 


77 


84 


Total length 


apical 


M 


7327 


5645 


4882 






MAD 


489 


637 


853 






IQR 


846 


1398 


1568 




basal 


M 


4439 


3398 


3664 






MAD 


518 


672 


1461 






IQR 


945 


1461 


1263 



Generated morphologies can be generated in isolation or in a forest setting. Basal 
and apical dendrites are treated separately. Presentation as in Table 3. Values 
shown for the morphologies from the Kawaguchi archive and the morphologies 
generated in "Isolation" and in the "Forest" setting (see main text). 

the neurons generated in isolation. In such cases, a neurite ter- 
minates once it leaves the designated space (Cuntz etal., 2010; 
Memelli et al., 2013). This might explain in part why the neurons 
generated in isolation and in a limited space show less variance 
(Figures 3H,I and 4E,F). However, since one of the future goals of 



this work is to generate full circuits, and because synapse occur- 
rence is proportional to structural overlap between axons and 
dendrites (Peters and Feldman, 1976), we cannot constrain the 
space and generate large ensembles of neurons simultaneously (as 
in the forest setting, Figures 3G and 4C). Therefore, future work 
will also focus on the design of proper rules for branching and 
termination. 

We can also start an argument by assuming that the exem- 
plar is not representative for all neurons. It has been demon- 
strated that reconstructed neurons contain a lot of biases related 
to reconstruction methods and selection by the experimenter 
(Horcholle-Bossavit etal., 2000; Kaspirzhny etal., 2002; Szilagyi 
and De Schutter, 2004; Steuber etal, 2004). For instance, the 
experimenter might select only "typical" neurons that are labeled 
well in the slice, which leads to a strong bias in the data. Also, neu- 
rons at the edge of a slice are more likely to be selected for technical 
reasons while precisely these neurons might be affected by the slice 
preparation in that neurites might be cut. Because these biases are 
not documented it is hard to make an estimate of their effect on 
the sample. As such, another option remains to explain the large 
variance in the generated data remains: the construction rules can 
be incomplete. Clearly the rules employed in this work are phe- 
nomenological and only crudely mimic morphological traits, so 
are incomplete. But assuming the rules are correct has interesting 
implications mainly because of the predictive power associated 
with a mechanistic model. Having a mechanistic explanation of 
neuron morphology has the advantage that morphological traits of 
various kinds can be predicted. For instance, age has an influence 
on morphologies and makes classifying neurons of varying age to 
correct classes nearly impossible (but see da Fontoura Costa et al., 
2002). With a mechanistic model, morphologies corresponding 
to a certain age could be generated and serve as ground truth. 
Similarly, to assess pathological cases, simulated knock-outs could 
be predicted. Predictions, in turn, could be used to validate the 
phenomenological construction rules: predict the outcome of a 
particular knock-out and compare the resultant traits in silico and 
in vitro. 

In conclusion, we designed, implemented and validated a new 
computational framework in accordance to a paradigm shift in 
the study of neuronal morphologies: away from studying mor- 
phologies in isolation to a study of neuronal morphologies as 
participants in their neuronal context. We demonstrated the 
potential of this new framework to study variation in neuronal 
morphology through a "generative" approach. Future research 
will focus on the generation and emergence of complete micro- 
circuits. 

REFERENCES 

Alvarez, F. J., Perason, J. C, Harrington, D., Dewey, D., Torbeck, L., and 
Fyffe, R. E. W. (1998). Boutons on alpha-motoneurons in the lumbar spinal 
cord of adult cats. /. Comp. Neurol. 393, 69-83. doi: 10.1002/(SICI)1096- 
9861(19980330)393:K69::AID-CNE7>3.0.CO;2-O 

Anwar, H., Riachi, I., Schurmann, F., and Markram, H. (2009). "An approach 
to capturing neuron morphological diversity," in Computational Neuroscience: 
Realistic Modeling for Experimentalists, ed. E. De Schutter (Cambridge: The MIT 
Press), 21 1-232. 

Ascoli, G. A., Donohue, D. E., and Halavi, M. (2007). NeuroMorpho.Org: a 
central resource for neuronal morphologies. /. Neurosci. 27, 9247-9251. doi: 
10.1523/JNEUROSCI.2055-07.2007 



Frontiers in Neuroanatomy 



www.frontiersin.org 



September 2014 | Volume 8 I Article 92 | 12 



Torben-Nielsen and De Schutter 



Context-aware modeling of neuronal morphologies 



Ascoli, G. A., Jeffrey, L. K, Ruggero, S., Slawomir, J. N., Stephen L. S, and 
Krichmar, G. L. (2001). Computer generation and quantitative morphome- 
tric analysis of virtual neurons. Anat. Embryol. (Berl.) 204, 283-301. doi: 
10.1007/s004290100201 

Banker, G. A., and Cowan, W. M. (1977). Rat hippocampal neurons in dispersed cell 
culture. Brain Res. 126, 397-425. doi: 10.1016/0006-8993(77)90594-7 

Cannon, R. C, Turner, D. A., Pyapali, G. K., and Wheal, H. V. (1998). An on-line 
archive of reconstructed hippocampal neurons. /. Neurosci. Methods 84, 49-54. 
doi: 10.1016/S0165-0270(98)00091-0 

Carim-Todd, L., Bath, K. G, Fulgenzi, G, Yanpallewar, S., Jing, D., Barrick, C. A., 
etal. (2009). Endogenous truncated TrkB.Tl receptor regulates neuronal com- 
plexity and TrkB kinase receptor function in vivo. /. Neurosci. 29, 678-685. doi: 
10.1523/JNEUROSCI.5060-08.2009 

Chen, J.-G., Rasin, M.-R., Kwan, K. Y., and Sestan, N. (2005). Zfp312 is required for 
subcortical axonal projections and dendritic morphology of deep-layer pyramidal 
neurons of the cerebral cortex. Proc. Natl. Acad. Sci. U.S.A. 102, 17792-17797. 
doi: 10.1073/pnas.0509032102 

Cove, J., Blinder, P., and Baranes, D. (2009). Contacts among non-sister dendritic 
branches at bifurcations shape neighboring dendrites and pattern their synaptic 
inputs. Brain Res. 1251, 30-41. doi: 10.1016/j.brainres.2008.1 1.028 

Cullheim, S., Fleshman, J. W., Glenn, L. L., and Burke, R. E. (1987). Membrane 
area and dendritic structure in type-identified triceps surae alpha motoneurons. 
/. Comp. Neurol. 255,68-81. doi: 10.1002/cne.902550106 

Cuntz, H., Forstner, F, Borst, A., and Hausser, M. (2010). One rule to grow them 
all: a general theory of neuronal branching and its practical application. PLoS 
Comput. Biol. 6:el000877. doi: 10.1371/journal.pcbi. 1000877 

da Fontoura Costa, L., Manoel, E. T., Faucereau, F, Chelly, J., van Pelt, J., and Ramak- 
ers, G. (2002). A shape analysis framework for neuromorphometry. Network 13, 
283-310. doi: 10.1088/0954-898X/13/3/303 

Dierssen, M., and Ramakers, G. J. (2006). Dendritic pathology in mental retarda- 
tiomfrom molecular genetics to neurobiology. Genes Brain Behav. 2, 48-60. doi: 
10. 1 1 1 II). 1601- 1 83X.2006.00224.X 

Dimitrova, S., Reissaus, A., and Tavosanis, G. (2008). Slit and Robo regulate dendrite 
branching and elongation of space-filling neurons in Drosophila. Dev. Biol. 324, 
18-30. doi: 10.1016/j.ydbio.2008.08.028 

Eberhard, J. P., Wanner, A., and Wittum, G. (2006). NeuGen: a tool for 
the generation of realistic morphology of cortical neurons and neural net- 
works in 3D. Neurocomputing 70, 327-342. doi: 10.1016/j.neucom.2006. 
01.028 

Frotscher, M. (1998). Cajal-Retzius cells, reelin, and the formation of layers. Curr. 

Opin. Neurobiol. 8, 570-575. doi: 10.1016/S0959-4388(98)80082-2 
Gao, F.-B. (2007). Molecular and cellular mechanisms of dendritic morphogenesis. 

Curr. Opin. Neurobiol. 17, 525-532. doi: 10.1016/j.conb.2007.08.004 
Graham, B. P., and van Ooyen, A. (2004). Transport limited effects in a model 

of dendritic branching. /. Theor. Biol. 230, 421-432. doi: 10.1016/j.jtbi.2004. 

06.007 

Grueber, W. B., Yang, C.-H., Ye, B., and Jan,Y.-N. (2005). The develop- 
ment of neuronal morphology in insects. Curr. Biol. 15, R730-R738. doi: 
10.1016/j.cub.2005.08.023 

Hentschel, H. G, and van Ooyen, A. (1999). Models of axon guidance and 
bundling during development. Proc. Biol. Sci. 266, 2231-2238. doi: 10.1098/rspb. 
1999.0913 

Hevner, R. E, Ray Daza, A. M., Rubenstein, J. L. R., Stunnenberg, H., Olavarria, J. 
E, and Englund, C. (2003). Beyond laminar fate: toward a molecular classifica- 
tion of cortical projection/pyramidal neurons. Dev. Neurosci. 25, 139-151. doi: 
10.1159/000072263 

Hill, S. L., Wang, Y, Riachi, I., Schurmann, E, and Markram, H. (2012). Statistical 
connectivity provides a sufficient foundation for specific functional connectivity 
in neocortical neural microcircuits. Proc. Natl. Acad. Sci. U.S.A. 109, E2885- 
E2894. doi: 10.1073/pnas.l202128109 

Hillman, D. (1979). "Neuronal shape parameters and substructures as a basis of 
neuronal form," in The Neurosciences, Fourth Study Program, ed. F. Schmitt 
(Cambridge, MA: The MIT Press), 477-498. 

Hintjens, P. (2013). ZeroMQ. O'Reilly Media. Available at: http://shop.oreilly.com/ 
product/0636920026136.do 

Hirai, Y, Morishima, M., Karube, E, and Kawaguchi, Y. (2012). Specialized corti- 
cal subnetworks differentially connect frontal cortex to parahippocampal areas. 
/. Neurosci. 32, 1898-1913. doi: 10.1523/JNEUROSCI.2810-11.2012 



Horcholle-Bossavit, G, Gogan, P., Ivanov, Y, Korogod, S., and Tyc-Dumont, S. 
(2000). The problem of the morphological noise in reconstructed dendritic 
asrborizations. /. Neurosci. Methods 95, 83-93. doi: 10.1016/S0165-0270(99) 
00159-4 

Itoh, T, Sobue, G, Yasuda, T, Kimata, K., Mitsuma, T, and Takahashi. A. 
(1993). Geometric response to nerve growth factor is preserved in aged rat sen- 
sory neurons: a single-neuron culture study. Neurobiol. Aging 14, 167-76. doi: 
10.1016/0197-4580(93)90093-0 

Jan, Y.-N., and Jan, L. Y. (2010). Branching out: mechanisms of dendritic 
arborization. Nat. Rev. Neurosci. 11, 316-328. doi: 10.1038/nrn2836 

Jan, Y.-N., Jan, L. Y, and Francisco, S. (2003). The control of dendrite development. 
Neuron AO, 229-242. doi: 10.1016/S0896-6273(03)00631-7 

Kaspirzhny, A. V., Gogan, P., Horcholle-Bossavit, G., and Tyc-Dumont, S. (2002). 
Neuronal morphology data bases: morphological noise and assesment of data 
quality. Network 13, 357-380. doi: 10.1088/0954-898X/13/3/307 

Kaufmann, W. E., and Moser, H. W. (2000). Dendritic anomalies in disor- 
ders associated with mental retardation. Cereb. Cortex 10, 981-991. doi: 
10.1093/cercor/10.10.981 

Koene, R., Tijms, B., van Hees, P., Postma, E, de Ridder, A., Ramakers, G, etal. 
(2009). NETMORPH: a framework for the stochastic generation of large scale 
neuronal networks with realistic neuron morphologies. Neuroinformatics 7, 195— 
210. doi: 10.1007/sl2021-009-9052-3 

Kriegstein, A., and Dichter, M. A. (1984). Neuron generation in dissociated cell 
cultures from fetal rat cerebral cortex. Brain Res. 295, 184-189. doi: 10.1016/0006- 
8993(84)90829-1 

Landgraf, M., and Evers, J. F. (2005). Control of dendritic diversity. Curr. Opin. Cell 

Biol. 17,690-696. doi: 10.1016/j.ceb.2005.09.005 
Lindsay, K. A., Maxwell, D. J., Rosenberg, J. R., and Tucker, G. (2007). A new 

approach to reconstruction models of dendritic branching patterns. Math. Biosci. 

205,271-296. doi: 10.1016/j.mbs.2006.08.005 
London, M., and Hausser, M. (2005). Dendritic computation. Annu. Rev. Neurosci. 

28, 503-532. doi: 10.1146/annurev.neuro.28.061604.135703 
Luczak, A. (2006). Spatial embedding of neuronal trees modeled by diffusive 

growth. /. Neurosci. Methods 157, 132-141. doi: 10.1016/j.jneumeth.2006. 

03.024 

Marin-Padilla, M. (1998). Cajal-retzius cells and the development of the neocortex. 

Trends Neurosci. 21, 64-71. doi: 10.1016/S0166-2236(97)01 164-8 
Marks, W. B., and Burke, R. E. (2007). Simulation of motoneuron morphology 

in three fimensions. i. building individual dendritic trees. /. Comp. Neurol. 503, 

685-700. 

McAllister, A. K. (2000). Cellular and molecular mechanisms of dendrite growth. 

Cereb. Cortex 10, 963-73. doi: 10.1093/cercor/10.10.963 
Memelli, H, Torben-Nielsen, B., and Kozloski, J. (2013). Self-referential forces 

are sufficient to explain different dendritic morphologies. Front. Neuroinform. 

7:1-12. doi: 10.3389/fninf.2013.00001 
Meyer, H. S., Wimmer, V. C, Hemberger, M., Bruno, R. M., de Kock, C. P., Frick, A., 

etal. (2010). Cell type-specific thalamic innervation in a column of rat vibrissal 

cortex. Cereb. Cortex 20, 2287-303. doi: 10.1093/cercor/bhq069 
Mortimer, D., Fothergill, T., Pujic, Z., Richards, L. J., and Goodhill. 

G. J. (2008). Growth cone chemotaxis. Trends Neurosci. 31, 90-98. doi: 

10.1016/j.tins.2007.11.008 
Oberlaender, M., de Kock, C. P. J., Bruno, R. M., Ramirez, A., Meyer, H. S., Dercksen, 

V. J., etal. (2012). Cell type-specific three-dimensional structure of thalamocor- 
tical circuits in a column of rat vibrissal cortex. Cereb. Cortex 22,2375-2391. doi: 

10.1093/cercor/bhr317 
Peters, A., and Feldman, M. L. (1976). The projection of the lateral geniculate 

nucleus to area 17 of the rat cerebral cortex. I. General description. /. Neurocytol. 

5, 63-84. doi: 10.1007/BF01176183 
Peters, A., and Payne, B. R. (1993). Numerical relationships between geniculocortical 

afferents and pyramidal cell modules in cat primary visual cortex. Cereb. Cortex 

3, 69-78. doi: 10.1093/cercor/3.1.69 

Samsonovich, A. V., and Ascoli, G. A. (2003). Statistical morphological analysis of 
hippocampal principal neurons indicates cell-specific repulsion of dendrites from 
their own cell. /. Neurosci. Res. 71, 173-187. doi: 10.1002/jnr.l0475 

Scott, E. K., and Luo, L. (2001). How do dendrites take their shape? Nat. Neurosci. 

4, 359-365. doi: 10.1038/86006 

Senft, S. L., and Ascoli, G. A. (1999). "Foundations and tools for neural modeling," in 
Proceedings of the International Work-Conference on Artificial and Natural Neural 



Frontiers in Neuroanatomy 



www.frontiersin.org 



September 2014 | Volume 8 I Article 92 | 13 



Torben-Nielsen and De Schutter 



Context-aware modeling of neuronal morphologies 



Networks, IWANN '99, Alicante, Vol. I, Lecture Notes in Computer Science 1606 
(Alicante: Springer), 25-33. 
Silver, R. A. (2010). Neuronal arithmetic. Nat. Rev. Neurosci. 11, 474-489. doi: 
10.1038/nrn2864 

Soltesz, I. (2005). Diversity in the Neuronal Machine: Order and Vari- 
ability in Interneuronal Microcircuits. New York, USA: Oxford University 
Press. 

Steuber, V., De Schutter, E., and Jaeger, D. (2004). Passive models of neurons in the 

deep cerebellar nuclei: the effect of reconstruction errors. Neurocomputing 58-60, 

563-568. doi: 10.1016/j.neucom.2004.01.095 
Stiefel, K. M., and Sejnowski, T. J. (2007). Mapping function onto neuronal 

morphology. /. Neurophysiol. 98, 513-26. doi: 10.1152/jn.00865 
Szilagyi, X, and De Schutter, E. (2004). Effects of variability in anatomical 

reconstruction techniques on models of synaptic integration by dendrites: a 

comparison of three internet archives. Eur. J. Neurosci. 19, 1257-1266. doi: 

10.1 1 1 l/j.l460-9568.2004.03222.x 
Tavosanis, G. (2012). Dendritic structural plasticity. Dev. Neurobiol. 72, 73-86. doi: 

10.1002/dneu.20951 

Torben-Nielsen, B., and Cuntz, H. (2014). "Introduction to dendrite morphology," 
in The Computing Dendrite: From Structure to Function, Vol. 11, Springer Series 
in Computational Neuroscience, eds H. Cuntz, W. Michiel, H. Remme, and B. 
Torben-Nielsen (New York, NY: Springer New York), 3-22. doi: 10.1007/978-1- 
4614-8094-5 

Torben-Nielsen, B., and Stiefel, K. M. (2010). An inverse approach for elucidating 
dendritic function. Front. Comput. Neurosci. 4:128. doi: 10.3389/fncom. 
2010.00128 

Torben-Nielsen, B., Vanderlooy, S., and Postma, E. (2008). Non-parametric algo- 
rithmic generation of neuronal morphologies. Neuroinformatics 6, 257-277. doi: 
10.1007/sl2021-008-9026-x 

van Ooyen, A., Carnell, A., de Ridder, S., Tarigan, B., Mansvelder, H. D., Bijma, F., 
etal. (2014). Independently outgrowing neurons and geometry-based synapse 
formation produce networks with realistic synaptic connectivity. PLoS ONE 
9:e85858. doi: 10.1371/journal.pone.0085858 



van Pelt, J., Carnell, A., de Ridder, S., Mansvelder, H. D., and van Ooyen. A. 

(2010). An algorithm for finding candidate synaptic sites in computer generated 

networks of neurons with realistic morphologies. Front. Comput. Neurosci. 4:148. 

doi: 10.3389/fncom.2010.00148 
van Pelt, J., and Schierwagen, A. (2004). Morphological analysis and modeling of 

neuronal dendrites. Math. Biosci. 188, 147-155. doi: 10.1016/j.mbs.2003.08.006 
van Pelt, J., and Uylings, H. B. M. (2002). Branching rates and growth functions 

in the outgrowth of dendritic branching patterns. Network 13, 261-281. doi: 

10.1088/0954-898X/13/3/302 
van Pelt, J., and Verwer, R. W. H. (1986). Topological properties of binary trees grown 

with order-dependent branching probabilities. Bell. Math. Biol. 48, 197-21 1. doi: 

10.1007/BF02460023 

Zubler, E, Hauri, A., Pfister, S., Bauer, R., Anderson, J. C, Whatley, 
A. M., etal. (2013). Simulating cortical development as a self construct- 
ing process: a novel multi-scale approach combining molecular and phys- 
ical aspects. PLoS Comput. Biol. 9:el003173. doi: 10.1371/journal.pcbi. 
1003173 

Conflict of Interest Statement: The authors declare that the research was conducted 
in the absence of any commercial or financial relationships that could be construed 
as a potential conflict of interest. 

Received: 07 June 2014; accepted: 20 August 2014; published online: 05 September 
2014. 

Citation: Torben-Nielsen B and De Schutter E (2014) Context-aware modeling of 
neuronal morphologies. Front. Neuroanat. 8:92. doi: 10.3389/fnana.2014.00092 
This article was submitted to the journal Frontiers in Neuroanatomy. 
Copyright © 2014 Torben-Nielsen and De Schutter. This is an open-access article 
distributed under the terms of the Creative Commons Attribution License (CC BY). 
The use, distribution or reproduction in other forums is permitted, provided the original 
author( s)or licensor are credited and that the original publication in this journal is cited, 
in accordance with accepted academic practice. No use, distribution or reproduction is 
permitted which does not comply with these terms. 



Frontiers in Neuroanatomy 



www.frontiersin.org 



September 2014 | Volume 8 I Article 92 | 14 



